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Motivated by recent studies reporting the formation of localized magnetic moments in doped 
graphene, we investigate the energetic cost for spin polarizing isolated impurities embedded in this 
material. When a well-known criterion for the formation of local magnetic moments in metals is 
applied to graphene we are able to predict the existence of magnetic moments in cases that are in 
clear contrast to previously reported Density Functional Theory (DFT) results. When generalized 
to periodically repeated impurities, a geometry so commonly used in most DFT-calculations, this 
criterion shows that the energy balance involved in such calculations contains unavoidable contribu- 
tions from the long-ranged pairwise magnetic interactions between all impurities. This proves the 
fundamental inadequacy of the DFT-assumption of independent unit cells in the case of magneti- 
cally doped low-dimensional graphene-based materials. We show that this can be circumvented if 
more than one impurity per unit cell is considered, in which case the DFT results agree perfectly 
well with the criterion-based predictions for the onset of localized magnetic moments in graphene. 
Furthermore, the existence of such a criterion determining whether or not a magnetic moment is 
likely to arise within graphene will be instrumental for predicting the ideal materials for future 
carbon-based spintronic applications. 

PACS numbers: 



The controlled introduction of impurities into an other- 
wise pristine solid is one effective way of tailoring the elec- 
tronic properties of a material. In systems such as thin 
films, nanowires and nanoparticles the effect of added im- 
purities is expected to be particularly pronounced due to 
the reduced dimensionality of the host materials. From a 
theoretical point of view, predictions of how the physical 
properties of a system are affected by dopants are cur- 
rently made based on state-of-the-art density functional 
theory calculations (DFT). These calculations often con- 
sider single impurities added to a unit cell with periodic 
boundary conditions, under the assumption that cells are 
sufficiently large and that impurities are not able to in- 
teract with their neighboring counterparts. There are, 
however, interactions that tend to become more long- 
ranged as the dimensionality is reduced, one of which 
is the interaction between magnetic impurities embed- 
ded in a metallic system. This raises the question of 
whether representing independent particles by single im- 
purity unit cells is a valid assumption in the case of mag- 
netic dopants in low-dimensional metallic structures. 

One example in which this assumption is frequently 
made is in the case of magnetic objects in graphene- 
related materials, of growing interest lately due to po- 
tential use in spintronic applications. As a matter of 
fact, a recent study has comprehensively investigated the 
magnetic properties of transition-metal atoms embedded 
in a graphene sheet [l| indicating complex magnetic be- 
havior as one moves across the periodic table. One re- 
markable finding in this pioneering survey is the absence 
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of a magnetic moment when one Fe atom substitution- 
ally replaces a single carbon atom of the graphene sheet. 
Of all transition metal atoms, it is somewhat surprising 
that such an iconic magnetic element like Fe seems un- 
able to develop a magnetic moment when immersed in 
graphene. Because the aforementioned DFT assumption 
of independent unit cells is used in the referred survey as 
well as in other studies of magnetic dopants in carbon- 
based structures ^ jjj , it is instructive to ask whether 
the intrinsic long-ranged interaction that arises between 
magnetic moments in low-dimensional metals might be 
responsible for interfering with some of the results re- 
cently reported. If so, this interference may spuriously 
suppress the formation of magnetic moments where they 
should actually exist. 

Lieb's theorem [S;] is often quoted to explain the mag- 
netic properties of graphene. It shows that a net mag- 
netization arises when there is an imbalance between the 
two non-equivalent sub-lattices composing the bipartite 
lattice of graphene. While this is a perfectly sound ex- 
planation for graphene flakes and ribbons, as well as 
for graphene sheets containing vacancies [6|, it is not 
directly applicable to substitutionally doped impurities 
since Lieb's theorem assumes a homogeneous electron- 
electron interaction throughout the system. In partic- 
ular, the considerably narrower d-band associated with 
transition-metal impurities makes the electronic interac- 
tion highly non-homogeneous, and another explanation 
for the origin of magnetic moments in doped graphene- 
related materials is required. The formation of a single 
local moment in a non-magnetic system has been gen- 
erally addressed by several authors, and a criterion for 
its existence has been previously derived in different con- 
texts ■ Here we sketch an alternative derivation of this 
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criterion, and generalize it to a pair of impurities in or- 
der to clarify the role played by the long range interaction 
between magnetic moments in low-dimensional systems. 
Furthermore, we show how this may affect first-principles 
calculations which artificially assume, for computational 
purposes only, that the system is translationally invari- 
ant. Although our focus is on doped graphene sheets, our 
conclusions result essentially from the hexagonal symme- 
try of the underlying lattice and are valid to graphene 
ribbons and flakes as well as nanotubes. 

We start by considering a single transition metal atom 
embedded in a non-magnetic host, which in our case 
is a pristine hexagonal lattice. We describe the elec- 
tronic structure of the system by a Hubbard-like Hamil- 
tonian H = Ho+ Hint, where Ho = 7,7 CjV<t 

represents the electronic kinetic energy plus a spin- 
independent local potential, and Hint is the electron- 
electron interaction term. The operator c|^^ creates an 
electron with spin a in atomic orbital /x on site i. We 



assume that Hint is an on-site interaction which takes 
place between electrons occupying the d orbitals of the 
transition metal impurity only, and is neglected else- 
where. In this case, the matrix elements of the spin- 
dependent part of the one-electron hamiltonian reduces 
to v'^j, — —^A^Sf_i^a, where represents the local ex- 
change splitting associated with orbital fj,, and a = ±1 
for 1 and J, spin, respectively. For simplicity we shall as- 
sume that the on-site effective exchange integrals U are 
the same for all d orbitals, whence it follows that A^ is 
/i-independent and equal to A = UM, where M is the 
local magnetic moment. 

To derive a local moment criterion we examine the sta- 
bility of the non-magnetic state described by Hq when 
Hint is activated considering relatively small values of 
the exchange splitting. The energy cost involved in the 
formation of a local magnetic moment at the impurity 
site is given by 



A£i = -ImTr / dE 

TT 



log 
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where GofiiE) is the single-particle Green function for 
an electron with energy E at the impurity site and Ep 
is the Fermi energy. The trace operator is over the or- 
bital degrees of freedom which in this case can be lim- 
ited to the five d-orbitals. The first term of the inte- 
grand in Eq.((T]) represents the reduction of the effective 
electron-electron interaction due to the appearance of a 
local spin imbalance at the impurity site, and the sec- 
ond term accounts for the increase of the correspond- 
ing electronic kinetic energy. The sign of A£i deter- 
mines whether or not the non-magnetic state is unsta- 
ble to a local magnetic moment formation. To derive 
a criterion for such an instability it is sufficient to ex- 
pand Eq.([T]) in powers of A to lowest order, which sim- 
plifies to A£i = {^UfiEp) + £(£;f)} (A/2)2, where 

({Ep) = ^ J^^dElmTi [Go.oiE)]'^ is the local suscepti- 
bility. The formation of a local magnetic moment at the 
impurity is then energetically favorable when A£i < 0, 
i.e., when 

KEf) > ^ . (2) 

This inequality sets the condition for the spontaneous for- 
mation of a single localized magnetic moment in a non- 
magnetic host. Written in terms of single-particle Green 
functions, it is model-independent and can be evaluated 
once the Hamiltonian is fully specified. The same crite- 
rion has been previously derived in other circumstances 
0, but here we have obtained it by total energy bal- 
ance considerations because it provides an easier way to 



analyze the effect of more impurities. 

Let us now imagine that a second transition-metal im- 
purity is added at site to. Similar steps may be taken 
to derive the following expression for the energy cost in- 
volved in the formation of local magnetic moments in the 
two impurities 

A2 /"^^ 

A£2 = 2 A£i + — / dElmTr [GqME) Gm.oiE)] . 

(3) 

The energy cost A£2 is not twice as large as A£i due 
to the interference between the two impurities. This 
is evident in the second term of Eq.([3]), which contains 
Green function propagators between sites and to. Most 
remarkably, this interference term that arises naturally 
when we evaluate the energy cost for the formation of 
two separate magnetic moments is precisely the same as 
the expression for the RKKY coupling between magnetic 
impurities embedded in a metallic system d, Q . Bearing 
in mind that this coupling is negative (positive) when the 
magnetic moments are parallel (antiparallel) and that it 
decays very slowly with the impurities separation in the 
case of low dimensional systems, this additional interfer- 
ence term may have striking consequences to the criterion 
presented above. 

Consider for instance a hypothetical impurity that 
meets the inequality of Eq.(l2]), that is, an impurity that 
possesses a magnetic moment when immersed in the 
graphene lattice. Suppose that we add a second impu- 
rity of the same hypothetical element with the imposed 
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constraint that both moments must be parallel to each 
other. The energy cost that was negative for a single 
impurity may become positive if the RKKY coupling 
favors an antiparallel alignment between the magnetic 
moments. In this case, may become positive even 
though Afi < 0. If this occurs, the two magnetic impu- 
rities whose moments are forced to remain parallel may 
adopt an altogether non-magnetic configuration rather 
than the most favorable antiparallel alignment. This is a 
clear indication that the artificial imposition of parallel 
alignment may introduce spurious effects as far as the 
determination of the true ground state configuration is 
concerned. 

DFT-based calculations that consider one single mag- 
netic impurity per periodically repeated unit cell implic- 
itly impose that their magnetic moments, should they 
exist, must be mutually parallel. Because of the peri- 
odic boundary conditions, the energy cost (per impu- 
rity) A£]^/N for inducing the spin splitting of iV equally 
spaced impurities becomes 

AF A? ^ f^" 

= ^^i + 4]V^ E J dElraTv [Go^njiE] Gn^^AE)] 

... (4) 

In this case the correction to the single-impurity con- 
tribution A£i, which once again is assumed to be neg- 
ative, is a sum of terms proportional to the pairwise 
RKKY interactions that may be positive and sufficiently 
large to reverse the sign of AE^/N. One could ar- 
gue that the RKKY interaction, being traditionally os- 
cillatory as a function of separation, will alternate be- 
tween negative and positive terms in the sum that ap- 
pears in Eq.(l4]), which will then average out and never 
be able to reverse the sign imposed by A£i. While 
this may be true in general, for graphene-based materi- 
als the underlying hexagonal atomic structure introduces 
a peculiar feature in the RKKY-like coupling that will 
seldom vanish the sum of Eq.(l4]). It has been shown 
that the magnetic coupling between impurities embed- 
ded in graphene-related materials obeys the following 
rule %^\: Ja,a/\Ja,a\ = Jb,b/\Jb,b\ = -Ja,b /\Ja,b\, 
where J a, A {Jb,b) represents the Heisenberg-like mag- 
netic coupling between impurities occupying the A {B) 
sub-lattice sites of the hexagonal lattice. In other words, 
in the case of graphene-related materials the sign of the 
coupling between impurities occupying equivalent sites 
is the same regardless of their separation. Therefore, if 
the magnetic coupling between like sites favors the anti- 
parallel alignment between the magnetic moments, all 
terms appearing in the summation of Eq.(|3]) will be pos- 
itive, which may converge to a sufficiently large value 
capable of overturning the satisfied criterion for single 
impurities. In addition, in the case of nanotubes, we 
have shown that the coupling magnitude tends to decay 
rather slowly as 1/D, where D is the separation between 
magnetic impurities Q. Such a slow decaying rate will 
turn the summation of Eq.Q into a non-convergent se- 
ries in the limit N oo, meaning that the correction 



will always surpass the magnitude of ASi . 

The striking implications of this mathematical analy- 
sis are that spurious nonmagnetic solutions may be ob- 
tained if existing magnetic moments are artificially con- 
strained to adopt a parallel alignment when they would 
spontaneously prefer to be antiparallel. This raises the 
question whether the recently reported absence of mag- 
netic moments for Fe in graphene could be one such case 
1]. It is worth stressing that more than correcting an 
inaccuracy, this would be a convincing indication of the 
inadequacy of the assumption, commonly used in DFT- 
calculations, of independent unit cells when dealing with 
magnetic dopants in carbon-based structures. One sim- 
ple way of testing if the moment suppression is the re- 
sult of the artificial constraint imposed by the periodic 
boundary conditions of the DFT scheme is to include 
more than one magnetic impurity per unit cell and al- 
low them to adopt both parallel and anti-parallel align- 
ments. In fact, in what follows we present DFT results 
for calculations comprising two impurities per unit cell 
and compare those with the results for a single impurity. 

Our DFT-calculations have been been made with the 
generalized gradient approximation [ll| for the exchange- 
correlation term. TrouUier-Martins pseudo-potentials 
[T^ and double-zeta polarization atomic orbitals [l^ 
have been used. The calculations were made with pe- 
riodic boudary conditions and supercells comprising 4x4 
and 7x7 graphene primitive cells. In these cases, one car- 
bon atom of the graphene lattice was substituted by a 
single Fe atom impurity. The distance between an im- 
purity and its image in the adjacent unit cell is 9.98 and 
17.47 A for the 4x4 and 7x7 supercells, respectively. Un- 
surprisingly, our results are very similar to those previ- 
ously reported We find that, when relaxed, the metal 
impurities are displaced outwards from the graphene sur- 
face by 1.14 A and that no magnetic moment is observed. 

The results are completely different, however, when 
two impurities per unit cell are considered. To maintain 
the same impurity separation as before, we duplicate the 
unit cells along one direction (4x8 and 7x14 primitive 
cells). With two Fe impurities per unit cell, we have the 
freedom to start these calculations with magnetic mo- 
ments in the anti-parallel configuration. Such an antifer- 
romagnetic alignment between the Fe moments is stable 
and energetically favorable, by 0.03 eV, when compared 
to the non-magnetic solution, which can be obtained by 
relaxing the spin-polarization. Remarkably, the substi- 
tutional Fe impurity does have a magnetic moment that 
is as large as 0.99 ^b- It is also instructive to calcu- 
late the energetics of the system in the configuration in 
which the Fe moments are parallel 14]. In this ferromag- 
netic case the total energy is considerably higher than 
the antiferromagnetic configuration. Table I shows the 
total energy values obtained for Fe as well as Mn impu- 
rities in the ferromagnetic (FM), spin-unpolarized (SU) 
and antiferromagnetic (AF) configurations. For Fe, the 
total energies in descending order are (FM, SU, AF). Be- 
cause the AF configuration is impossible to obtain with 
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the single-impurity unit cell, the system adopts the next 
possible configuration, which shows no spin polarization. 
Alternatively, one can understand this in terms of Eq. Q , 
which means that the magnetic-coupling correction that 
arises due to the periodic boundary conditions is able 
to revert the sign imposed by Af i leading to a spurious 
suppression of the existing Fe magnetic moments. 





Fe 


Mn 


FM 


0.14 


0.04 


SU 


0.03 


1.80 



TABLE I: Total energies, in eV, of the FM and SU configu- 
rations, for Fe and Mn impurities in graphene. All quantities 
are expressed relatively to the total energy of the AF con- 
figuration, which is the most energetically favorable for both 
impurities. 

Also shown in Table I are the values for Mn impuri- 
ties, for which the total energies in descending order are 
(SU, FM, AF). Another interesting result, not considered 
in reference [l|, is that the AF is the most energetically 
favorable configuration. Once again, this is easily under- 
stood by the single-impurity unit cell constraint that is 
unable to account for the AF alignment of the magnetic 
moments. In this case, the next possible configuration 
is the FM alignment, which again explains the results of 
reference In terms of Eq.Q, the magnetic-coupling 
correction for Mn impurities are not sufficient to overturn 
the sign determined by i , which means that A£n /N is 
still negative justifying the splitting of the spin-polarized 
bands into a FM configuration. 

This is an unmistakable proof of the potential prob- 
lems that may arise when dealing with magnetic impu- 
rities in graphene-based structures through the standard 
DFT scheme of single-impurity unit cells. Furthermore, 
it is also a clear indication of the relevance of the coupling 
between mag netic impurities across graphene-related ma- 
terials H, [la El m, O m, El ■ When this coupling is 
positive and sufficiently large to reverse the sign imposed 
by A£i , the artificial constraints imposed by the periodic 
boundary conditions spuriously suppresses the magnetic 
moment that would spontaneously exist in isolation. In 
DFT-calculations of doped graphene-related materials, it 
is therefore of paramount importance to consider more 
than a single impurity per unit cell and study the ener- 
getics of all possible configurations, namely, FM, SU and 
AF. 

Finally, regarding the condition for the formation of a 



localized magnetic moment expressed by the inequality of 
Eq.®, we can test its predictive power by applying it to 
the cases considered here. Written in terms of single- 
particle Green functions, the susceptibility i{Ep) can 
be further simplified in the case of small spin-splittings 
(A -C 1) to 1{Ef) ~ pq{Ef), where po{Ep) is the spin- 
unpolarized local density of states (LDOS) at the impu- 
rity site evaluated at the Fermi level Ep- The value of 
U ~ leV, being primarily an atomic property, is fairly 
constant for all transition metal elements [2l|. Therefore, 
we can use the LDOS obtained at our SU calculations and 
test whether the inequality of Eq. Q is satisfied. Reassur- 
ingly, Upo{Ep) > 1 for both Fe and Mn, indicating that 
both elements favor the formation of a magnetic moment 
when embedded within graphene. Further tests were car- 
ried out with Ni impurities. In this case, the low value 
found for pq(Ep) does not meet our criterion, suggesting 
that Ni atoms within graphene will not develop a mag- 
netic moment. In fact, this is what we have found in our 
DFT-calculations with two Ni impurities per unit cell, 
which also agrees with previously reported results [11, [2] ■ 
Such a good agreement with the predictions based on 
Eq.® indicates that SU calculations, which are consid- 
erably less time consuming than spin-polarized ones, can 
be carried out first to test whether a localized magnetic 
moment is likely to arise. If so, further spin-polarized cal- 
culations arc required in which all possible configurations 
must be then considered. 

In summary, we have shown that the use of single- 
impurity-doped unit cells in DFT-based calculations 
is highly inappropriate to describe magnetically doped 
graphene and that it may lead to fundamentally erro- 
neous results. This is a consequence of the inherently 
long range nature of the magnetic interaction between 
impurities, which makes the hypothesis of independent 
unit cells in such systems invalid. Graphene being a ma- 
terial of increasing popularity and DFT-calculations be- 
ing the most widely used tool for studying the physical 
properties of materials, it is of paramount importance to 
account for this interaction when describing magnetically 
doped graphene and related structures. Furthermore, we 
present a mathematically transparent criterion for the 
formation of magnetic moments in graphene, something 
that has only been previously attempted on an ad-hoc 
basis The existence of a simple criterion that can 
tell whether or not a magnetic moment will arise when 
impurities are introduced to graphene-related materials 
is essential to predict which of these structures may be 
useful for spintronic applications. 
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